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An empirical extreme value distribution 
with tower and upper bounds proposed 
by the author is applied to represent 
probability distribution models for max- 
imum load intensities of the earthquake 
ground motion, the wind speed, and 
the live load in supermarkets. One of 
the difficulties in the estimation of the 
parameters is determining the upper 
bound value. Nevertheless application 
of Ihe proposed distribution to the an- 
nual maximum earthquake ground mo- 
tion results in considerable improve- 
ments over other models. Possible im- 



provements to the annual maximum 
wind speed model are discussed. The 
proposed distribution is also a good 
candidate for the live load extremes. 
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1. Introduction 

Probabilistic load models are utilized for limit 
state design procedures and safety assessments of 
structures. Since lifetime maximum loads have to 
be applied to these analyses, appropriate probabil- 
ity distributions are needed to represent load in- 
tensity models. The Gumbel distribution (Type I 
extreme value) and the Frechet distribution (Type 
II extreme value) are often used for such purposes. 

When the coefficient of variation (cov) is not 
large, discrepancies of the upper tails may not be 
very serious. However for load intensities with 
fairly large cov such as earthquake ground motions, 
existing extreme value distributions do not provide 
good fits to statistical data. Then an empirical ex- 
treme value distribution with both upper and lower 
bounds proposed by the author [1] is a good alter- 
native to improve probabilistic load models. 

Statistical data were prepared for the annual 
maximum earthquake ground motion, the annual 



maximum wind speed, and live load extremes in- 
crowding situations. The significance of the pro- 
posed distribution is discussed in terms of lifetime 
maximum statistics. 

2. Proposed Extreme Value Distribution 

Three types of extreme value distributions are 
commonly used for engineering purposes. Cumula- 
tive distribution functions-are written as follows [2], 

Fi = exp[-exp {-a (*-£>)}] - » <jc < » (1) 
F„(x)=exp[-(^) T ] *<x<oo (2) 



**>««*[-e^)'] 



- <» <jc <w (3) 
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where a , b, c, e, y, w, v are parameters which char- 
acterize the form of distribution. Distributions ex- 
pressed in Eqs. (1), (2), and (3) are the Gumbel, 
Frechet and Weibull distributions respectively. In 
Eq. (1), the random variable x could theoretically 
vary between - o° and + oo, while in Eq. (2) the 
lower bound value, e, and in Eq. (3) the upper 
"bound value, w, exist. When natural phenomena are 
considered, it seems reasonable that the physical 
quantity has a positive value with an upper bound 
limit. On consideration the formula of Eqs. (2) and 
(3), the following empirical extreme value distribu- 
tion has been proposed fl] 



W = exp[-k^} T ], 



e<x <«r 



(4) 



where w and e are upper and lower bound values, 
respectively, and u and y are scale and shape 
parameters, respectively. It can be seen that whenx 
approaches the lower bound, e, Eq. (4) approaches 
Eq. (2) with c=(w - e)/u and when x approaches 
the upper bound, w, Eq. (4) approaches Eq. (3) 
with v = w -u(w — e). 

In order to demonstrate the form of the proposed 
distribution, the effect of parameter u with w = 10 
and y = 1.0 for Eq. (4) is shown in Fig. 1 on Gumbel 
probability paper, with x on the ordinate and re- 
duced variate y on the abscissa. In a similar way the 
effect of parameter y with w = 10 and u =exp(4/>') 
for Eq. (4) is demonstrated in Fig. 2. 



10- 



x 

a> 

CO 

> 

e 

3 

E 

1 

E 



■■ 1 


i 1 i i i i r — i — i — i— 


i i 


//^/^/^ 


jf - 


- 


**■// if / / / 


- 




MSN J 

/ / 7 ™/#/ 

/ / / / II/sl 
/ / / / / / 7 

tf^t^-tfrf^t-^i i i i i 


i i 



-202468 
reduced variate, y 



10 



> 
E 

E 

s 




2 4 6 

reduced variate, y 



Fig. 2. Effect of parameter y on proposed distribution with 
h> = 10, u =exp(4/y). (Gumbel distribution probability paper.) 



3. Annual Maximum Earthquake 
Ground Motion Model 

The seismic hazard estimation has often been 
based on earthquake occurrence models assuming a 
Gutenberg-Richter relationship. An alternative es- 
timation is possible when sufficient number of 
earthquake records are available to acquire annual 
maximum earthquake ground motion data at a site. 
Such an approach is rather common in Japan begin- 
ning with Kawasumi's work in 1951 [3]. A recent at- 
tempt was made by applying the proposed 
distribution of Eq. (4) [4]. Some modifications were 
introduced in this study. Earthquake data for the 
last 400 years were utilized according to Usami's 
catalogue [5]. Kanai's attenuation law was chosen 
as a representative relationship between the 
bedrock velocity, V, and the magnitude, M, with the 
hypocentral distance,*, expressed in the following 
formula [6] 



V ■=■ irjMi-w-ci.Mi+i^iogi-to.oi+Ai?) 



(5) 



Fig. 1. Effect of parameter u on proposed distribution with 
h> = 10, y=*1.0. (Gumbel distribution probability paper.) 



where a focal depth of 30 km is uniformly assumed 
to calculate*. 

The annual maxima of bedrock velocity calcu- 
lated according to Eq. (5) are plotted on Gumbel 
probability paper for four sites, i.e., Sendai, Tokyo, 
Osaka, and Fukuoka in Fig. 3. The 50 largest data 
from annual maxima in the 400 year period were 
used, since minor earthquake motions were consid- 
ered to be missing in historical records or docu- 
ments and so should be eliminated from the 
analysis. 
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Fig. 3. Extreme value fitting to annual maximum bedrock velocity of earthquake motion in Japan, where the solid and dashed lines 
indicate the proposed and the Frcchct distribution, respectively. 



The upper bound value could be assumed from 
tektonic findings on the fault activity [4], however 
the upper bound magnitude seems to provide a 
fairly rough estimate for a particular site as a varia- 
tion of magnitude by 0.1 causes only 15% change in 
the estimation of V. 



The proposed distribution fitted to plotted data 
by the least squares method is shown where the 
value of w in Eq. (4) was chosen by engineering 
judgment as w =5.0, 12.0, 10.0 and 3.0 for Sendai, 
Tokyo, Osaka and Fukuoka, respectively. The 
Frechet distribution fitted by the same method is 
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also shown in a dashed line for comparison. The ex- 
istence of saturation tendency indicates that the 
representation by Eq. (4) is better. 

The difference in the results between the pro- 
posed distribution and the Frechet distribution can 
be summarized in Table 1. The error is estimated in 
terms of the normalized square root of sum of 
squares error as, 



£•-■■^2 <x, -*)»/£ 2*. 



(6) 



where x, is the ith annual maximum data and*, is the 
corresponding value estimated from the distribu- 
tion model, and n = 50 for earthquake models. 

Although the extimated mean of 50 year maxi- 
mum values based on the Frechet distribution is 
similar to that based on the proposed distribution, 
estimated cov values for the 50 year maxima based 
on the Frechet distribution, which are the same as 
those for annual maxima, are considerably greater 
than those yielded by the proposed distribution. 
Significant reduction in the error estimate also indi- 
cates the appropriateness of the proposed distribu- 
tion. The cumulative distribution of the 50 year 
maximum was obtained as the 50th power of the cu- 
mulative distribution function of the annual maxi- 
mum, i.e., it was assumed that the annual maxima 
are mutually independent. The mean and cov of the 
50 year maximum were calculated numerically for 
the proposed distribution as the closed form rela- 
tionship between the mean and variance and the 
parameters in Eq. (4) is not obtainable. 

4. Annual Maximum Wind Speed 

The Gumbel distribution is often used to repre- 
sent the annual maximum wind speed distribution. 



The possibility of improved representation by Eq. 
(4) is examined for sites where some saturation ten- 
dencies are observed, i.e., Aomori, Akita, Nagoya 
and Kagoshima. 

The wind speed data were corrected by taking 
into account changes of measurement height and 
the change of the terrain roughness in the period 
between 1960 and 1970 [7]. Measured data at 
meteorological agency stations in the period be- 
tween 1929 and 1991 were utilized. Plotted data and 
distribution curves fitted to the plots, as was done 
for the earthquake cases, are shown in Fig. 4, with 
dashed lines representing the Gumbel distribution. 
The upper bound w = 35 (m/s) was used for Aomori 
and Akita, while w =40 was used for Nagoya and 
Kagoshima in Eq. (4). Although the difference be- 
tween the two types of distributions is not as signif- 
icant as in the case of earthquakes, error estimates 
are improved except for Nagoya, where the fitting is 
rather poor in comparison with other cases as seen 
in Table 2. The use of a nonzero lower bound value, 
e.g., e «= 10, could improve the fitting for the case of 
Nagoya. However, this was avoided. Two different 
major factors, such as the occurrence of typhoons 
and monsoons, could be the reason for the concave 
shape of plots on the Gumbel probability paper. 

Estimated mean and cov values for the annual 
maxima and the 50 year maxima are also listed in 
Table 2. The reduction in estimation of cov of the 
50 year maximum for the proposed distribution can 
be pointed out. When the existence of an upper 
bound for the extreme value distribution of wind 
speed is accepted, such a reduction could result in 
a smaller load factor in the probability-based 
design. 



Tabic I. Statistics of maximum earthquake ground motion 



Site 


B" 


Proposed 
u 


distribution, Eq 
7 E 


(4) 
Annual 


max 


50 year 


max 


c 


Frechet dist 

y E 


ribution, 
Annual 


Eq.(2) 
max 50 


year max 




Mean 


cov 


Mean 


cov 


Mean 


cov 


Mean 


Scndai 


5.0 


28.7 


1.30 


0.11 


0.32 


1.56 


2.52 


0.33 


0.455 


2.88 


0.23 


0.63 


0.74 


2.44 


Tokyo 


12.0 


57.6 


1.15 


0.125 


0.61 


1.73 


5.46 


0.48 


0.660 


2.34 


0.25 


1.03 


1J0 


5.49 


Osaka 


10.0 


90.9 


1.18 


0.21 


0.34 


1.99 


3.48 


0.60 


0.303 


2.06 


0.37 


0.52 


3.29 


3.50 


Fukuoka 


3.0 


45.8 


1.19 


0.10 


0.18 


1.58 


1.43 


0.44 


0.199 


2.51 


0.24 


0J0 


1.03 


1.41 
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Fig. 4. Extreme value fitting to annual maximum wind speed in Japan, where the solid and dashed lines indicate the proposed and the 
Gumbcl distribution, respectively. 
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Table 2. 


Statistics of maximum 


wind speed 
























Site 


w 


Proposed distribution, 
u y E 


Eq. (4) 
Annual max 


50 year 


max 


a 


Gumbel d 
b E 


sirtbution, Eq. (1) 
Annual max 50 ycai 


max 




Mean 


COV 


Mean 


cov 


Mean 


cov 


Mean 


cov 


Amori 


35 


1.09 


3.00 


0.02 


18.4 


0.19 


27.9 


0.08 


037 


16.8 


0.02 


18.4 


0.19 


29.1 


0.12 


Akita 


35 


0.70 


2.98 


0.02 


22.0 


0.14 


30.1 


0.05 


0.41 


20.6 


0.03 


22.0 


0.14 


31.7 


0.10 


Nagoya 


40 


1.63 


2.79 


o.os 


17.3 


0.25 


29.8 


0.10 


0.29 


15.4 


0.08 


17.4 


0.25 


30.9 


0.14 


Kagoshima 40 


U9 


2.05 


0.02 


19.4 


0.29 


34.0 


0.08 


0.23 


16.9 


0.03 


19.4 


0.29 


36.8 


0,15 



5. Simulated Extreme Live Load 

Extraordinary live loads for the ultimate limit 
state design may be estimated based on computer 
simulations according to a scenario for extraordi- 
nary situations. One reported possibility is to model 
crowd-gathering situations in supermarkets [8]. 
Room plans with specified shapes and weights of 
racks with goods and other furnishings were used 
according to surveyed data, and crowd-gathering 
situations at one corner of each room were simu- 
lated for typical cases such as, 

(a) from 0.3 person/m 2 to 5.0 person /m 2 

(b) from 1.0 person/m 2 to 10.0 person /m 2 

Personnel loads were distributed in the area with 
no racks or furniture, and 700 Wwas postulated as 
the weight of a person. 



The equivalent uniformly distributed loads 
(EUDL) were calculated for slab end bending mo- 
ments in a shorter span and for girder end bending 
moments. The detailed procedure is described else- 
where [8]. 

The plotted data with distribution curves of Eqs. 
(1) and (4) are shown in Fig. 5 in a similar manner 
to Figs. 3 and 4. Dashed lines represent the Gumbel 
distribution of Eq. (1). The upper half of the data 
were used to obtain parameters of distributions by 
the least squares method. The upper bound values 
of 4 kPa and 7 kPa were used for cases (a) and (b), 
respectively. 

Results are summarized in Table 3. Saturation 
tendency is clear for case (b) where the error is sig- 
nificantly reduced by the proposed distribution. 
Assumed personnel load intensities in the scenario 
are somehow arbitrary and a further survey could 
improve models introduced herein. Nevertheless 
the usefulness of the proposed extreme value distri- 
bution with upper bound can be recognized. 



Table 3. Satisfies of extreme EUDL for supermarket 



Case 




Proposed distribution 


Eq.(4) 






Gumbel distribution, 


Eq.(l) 


w 


u 


y 


E 


Mean 


cov 


a 


b 


E 


Slab (a) 


4000 


3.65 


1.61 


0.04 


1674 


0.36 


000182 


827 


0.07 


Slab (b) 


7000 


13.6 


0.80 


0.11 


2476 


0.68 


0.00067 


193 


0.24 


Girder (a) 


4000 


4.22 


1.96 


0.03 


1390 


0.36 


0.00216 


676 


0.04 


Girder (b) 


6000 


14.3 


0.81 


0.10 


2075 


0.71 


0.00076 


64.7 


0.22 
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Fig. 5. Extreme value fitting to simulated extraordinary live loads due to crowding concentration in supermarkets, where the solid and 
dashed lines indicate the proposed and the Gumbcl distribution, respectively. 
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6. Concluding Remarks 

An empirical extreme value distribution with 
both upper and lower bounds was reviewed, The 
usefulness and improved fit to extreme load intensi- 
ties available, such as the annual maximum earth- 
quake ground motion, the annual maximum wind 
speed and the extreme EUDL due to crowding situ- 
ations in supermarkets, were demonstrated. Use of 
the simpler commonly used Gumbel or Frechet dis- 
tributions could cause some significant overestima- 
tion in the coefficient of variation for lifetime 
maximum loads. 
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